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Q ■ Abstract 

^ I We study the problem of discrete geometric packing. Here, given weighted regions 

0^ ' (say in the plane) and points (with capacities), one has to pick a maximum weight 

subset of the regions such that no point is covered more than its capacity. We provide a 

general framework and an algorithm for approximating the optimal solution for packing 
V^ ' in hypergraphs arising out of such geometric settings. Using this framework we get 

\^ • a flotilla of results on this problem (and also on its dual, where one wants to pick a 

t/3 ! maximum weight subset of the points when the regions have capacities). For example, 

for the case of fat triangles of similar size, we show an 0(l)-approximation and prove 

that no PTAS is possible. 
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^ ■ 1 Introduction 

• ' Motivation and examples. Consider the problem of obnoxious facility location 

(^ ■ |Tam91l[Cap99| ; that is, you have to place several facilities, but these facilities are undesired 

(i.e., obnoxious). Facilities of this type include nuclear reactors, wind farms, airports, power 

plants, factories, prisons, universities, etc. Facilities can also be semi-desirable - a customer 

might want to have supermarkets close to their home, but they do not want to have too 

/\ ' many of them close by as they increase traffic, noise, etc. One natural way to model this 

c^ ' geometrically is to associate each obnoxious facility with its region of undesirability. We 

also have customers (modeled as points), and each customer has a threshold of how many 

obnoxious facilities it is willing to accept covering it. Different customers may have different 

thresholds, for example because more affluent people have stronger political power and it is 

harder to place obnoxious facilities near their homes. 
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Naturally, if you allow only a single region to cover each customer, then this is a classical 
packing problem, and much work has been done on packing disks/balls |SMC"'"07] . However, 
there are many cases where allowing limited interaction between the packed regions is allowed 
(after all, these facilities are required for modern existence). As a concrete example of this 
type of problem, consider the placement of radio stations/cellphone towers. While airports 
allow only very limited levels of interference!, higher levels of such interference is acceptable 
in residential neighborhoods. However, at a certain point there is going to be resistance to 
placing more wireless towers in residential areas, as these towers are viewed as causing cancer 
(this fear might be baseless, but it does not change the political reality of the difficulty of 
placing such towers). On the other hand, there is little resistance to placing such towers 
along highways in sparsely populated areas. 

In this paper, we are interested in the modeling of such problems and in the computation 
of an efficient approximation to the optimal solution of such problems. 

Modeling. 

As hinted by the above, perhaps the most natural way to model this problem is as a gener- 
alization of the well known independent set problem. 

Independent set is a fundamental discrete optimization problem. Unfortunately, it is 
not only computationally hard, but it is even hard to approximate to within a factor of 
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for any constant e |Has99] (under the assumption that NP 7^ P). Surprisingly, the 
problem is considerably easier in some geometric settings. For example, there is a PTASO 
|Cha03t IEJS05] for the following problem: Given a set of unit disks in the plane, find a 
maximum cardinality subset of the disks whose interiors are disjoint. Furthermore, a simple 
local search algorithm yields the desired approximation: For any e > 0, the local search 
algorithm that tries to swap subsets of size 0(1/ e"^) yields a (1 — £:)-approximation in rP^^''^ > 
time pTOQllCHTT] . 

The discrete independent set problem. In this paper, we consider packing problems 
in geometric settings that are natural extensions of the geometric independent set problem 
described above. As a starting point, motivated by practical applications, we consider the 
discrete version of the geometric independent set problem in which, in addition to a set of 
weighted regions, we are given a set of points, and the goal is to select a maximum weight 
subset of the regions so that each point is contained in at most one of the selected regions. 
We refer to this problem as the discrete independent set problem. Chan and Har-Peled 
[CHll] studied this discrete variant and proved that one can get a good approximation if 
the union complexity of the regions is small. 



^See http : //tinyurl . com/7td67v3 for a story of an airport closing down because of radio interference. 
^Polynomial time approximation scheme. 





Note that the discrete independent set problem captures the 
continuous version of the independent set problem, since we can 
place a point in each face of the induced arrangement of the given 
regions. In fact, the discrete version is considerably harder (in 
some cases) than the continuous variant. The difficulty lies in 
that several regions forming a valid solution to an instance of a 
discrete independent set problem may contain a common point 
that is not part of the set of points given as input; the figure on 
the right shows an example in which the middle point, marked 
as a square, is covered twice by the given valid solution. 

To illustrate the difference in difficulty, consider the case when the 
input consists of a set S of segments (in general position) with their 
endpoints on a circle, such that every pair of segments intersect. Clearly, 
in the continuous version, the maximum independent set of segments is 
a single segment. However, in this case, the discrete version captures 
the graph independent set problem. More precisely, we can encode any 
instance of independent set (i.e., a graph G = {V,E)) as an instance 
of this problem as follows. Every vertex t; G l^ is mapped to a segment s„ of S, and every 
edge uv G E, is mapped to the point Su n s^, (which is added to a set of points P). Clearly, 
an independent set of segments of S (in relation to the point set P) corresponds to an 
independent set in G. That is, the geometric discrete version is sometimes as hard as the 
graph independent set problem. For example, the figure on the right depicts the resulting 
instance encoding independent set for K^^. 

The packing problem. In this paper, we are interested in the natural extension of the 
discrete independent set problem to the case where every point has a capacity and might be 
covered several times (but not exceeding its capacity). The resulting problem has a fiavor 
of a packing problem, and is defined formally as follows. 

Problem 1.1 (PackRegions.) Given a set T) of regions and a set P of points such that 
each region r has a weight w{r) and each point p has a capacity #(p), find a maximum weight 
subset X of the regions such that, for each point p, the number of regions in X that contain 
p is at most its capacity #(p). 

We emphasize that different points might have different capacities, which makes the 
problem considerably more challenging to solve than the unit capacities case (i.e., the discrete 
independent set problem). We also consider the following dual problem in which the points 
have weights and the regions have capacities. 

Problem 1.2 (PackPoints.) Given a set D of regions and a set P of points such that each 
region r has a capacity #(r) and each point p has a weight w{p), find a maximum weight 
subset X of the points such that each region r contains at most #(r) points of X. 



Hypergraph framework. These two problems can be stated in a unified way in the 
language of hypergraphqj. Given an instance of PackRegions, we construct a hypergraph 
as follows: Each weighted region is a vertex, and all the regions containing a given point of 
capacity k become a hyperedge (consisting of these regions) of capacity k. A similar reduction 
works for PackPoints, where the given weighted points are the vertices, and each region 
of capacity k becomes a hyperedge of capacity k consisting of all of the points contained in 
this region. Therefore the previous two problems are special cases of the following problem. 

Problem 1.3 (HGraphPacking.) Given a hypergraph G = (V, E) with a weight function 
w{-) on the vertices and a capacity function #(■) on the hyperedges, find a maximum weight 
subset X C V, such that \/f eE we have \X f] f\ < i^{f). 

We will be interested primarily in hypergraphs with certain hereditary properties. A 
hypergraph property is hereditary if the sub- hypergraph induced by any subset of the vertices 
has the property; an example of a hereditary property of hypergraphs is having bounded 
VC dimension. Roughly, we are interested in hypergraphs having the bounded growth 
property: For any induced sub- hypergraph on t vertices the number of its hyperedges that 
contain exactly k vertices is near linear in t and its dependency on k is bounded by 2'^^''\ see 
Definition 12.11 Such hypergraphs arise naturally when considering points and "nice" regions 
in the plane. 

Our results. 

• Main result. Our main result is an algorithm that provides a good approximation for 
HGraphPacking as a function of the growth of the hypergraph, see Theorem 13.111 Our 
result can be viewed as an extension of the work of Chan and Har-Peled |CH11] to these 
considerably more general and intricate settings. 

• Regions with low union complexity. In Section |U we apply our main result to regions that 
have low union complexity, and we get the following results: 

(A) If the union complexity of n regions is 0{nu{n)) then we get an O (u(?7,)^/'^) -approximation 
for PackRegions, where u is the minimum capacity of any point in the given in- 
stance. (That is, the problem becomes easier as the minimum capacity increases.) For 
the case where all the capacities are one, this is the discrete independent set problem, 
and our algorithm specializes to the algorithm of Chan and Har-Peled |CHllj , which 
gives an 0(M(?7,))-approximation. 

(B) More specifically, we get a constant factor approximation for PackRegions if the 
union complexity of the regions is linear. This holds for (i) fat-triangles of similar 
size, (ii) unit axis-parallel cubes in 3d, and (iii) pseudo-disks. See Corollary 14.31 

(C) Similarly, since the union complexity of fat triangles in the plane is 0{nlog* n) 

[EASlll lAdBESlTb] . we get an Oi{\og* n) '^\ approximation for such instances of 
PackRegions. 



^A hypergraph G is a pair (V, E), where V is a set of vertices and E is a collection of subsets of V which 
are called hyperedges . 



Bi-criteria approximation. Our main result also implies a bi-criteria approximation algo- 
rithm. That is, we can improve the quality of the solution, at the cost of potentially 
violating low capacity regions. Formally, if the input instance G = (V, E) of HGraph- 
Packing has at most Fk{t) = 2'-''^^^F{t) edges of size k when restricted to any subset 
of t vertices, then for any integer > 1, our algorithm yields an (0((-F(?7,)/n)^/'^) ,</>)- 
approximation to the given instance G of HGraphPacking. Specifically, the value of the 
generated solution X is at least r2(opt/(F(n)/n)^/''^), where opt is the value of the optimal 
solution, and for every hyperedge / G E, we have |/ fl X| < max(</), #(/)). 
As an example, for any set of n regions in the plane such that the boundaries of any 
pair of them intersects 0(1) times, the above implies that one can get an (0{v}^'l'^ , 0)- 
approximation for PackRegions. 

Axis-parallel boxes. The union complexity of axis-parallel rectangles can be as high as 
quadratic, and therefore we cannot immediately apply our main result to get a good 
approximation. Instead, we decompose the union of axis-parallel rectangles into regions 
of low union complexity, and this decomposition together with our main result gives us 
an O(logn) approximation for instances of PackRegions in which the regions are axis- 
parallel rectangles in the plane(see Lemma [4.9p . 

A more involved analysis also applies to the three dimensional case, where we get an 
O(log^n) approximation for PackRegions for axis parallel boxes (see Lemma [4. lip . 

Dual problem. We show in Section W2\ that, by standard lifting techniques, we can apply 
our result for PackRegions, where the regions are disks, to the dual problem of Pack- 
PointsInDisks. However, for other regions, the dual problem PackPoints seems to be 
more challenging. Specifically, this is true for the case of axis-parallel rectangles. For this 
case, we first provide a constant factor approximation for skyline instances of the problem; 
a skyline instance is a set of rectangles that lie on the x-axis. Interestingly, if the set of 
rectangles are defined in relation to a set of points (and each rectangles contains only a few 
points), then one can define a near- linear (in the number of points) sized set of rectangles 
such that each original rectangle is the union of two new rectangles. Combining this with 
the skyline result and a sparsifying technique, we get an(0(logn), 2)-approximation; that 
is, every rectangle b contains at most max(2, #(b)) points of the solution constructed, and 
the total weight of the solution is f2(opt/ log n) (see Theorem 14. 21 p . (Note that, by applying 
our general framework directly to this setting, we only get an {0{n^/'^^ , 0) -approximation, 
for any integer > 0.) 

Packing points into fat triangles. We provide a polylog bi-criteria approximation for the 
problem of packing points into fat triangles. This requires proving that one can compute, 
for a given point set, a small number of canonical subsets, such that the point set covered 
by any fat-triangle (if the set is sufficiently small), is the union of constant number of these 
canonical subsets. Proving this requires non-trivial modifications of the result of Aronov 
et al. |AES10j . In addition, we show that a measure defined over a fat triangle can be 
covered by a few fat triangles, each one of them containing only a constant fraction of 
the original measure. We believe these two results are of independent interest. Plugging 
this into the machinery, previously developed for axis parallel rectangles, yields the new 
approximation algorithm. See Section |5] for details. 

PTAS for disks and planes. We adapt the techniques of Mustafa and Ray |MR10j in order 



to get a PTAS for instances consisting of unweighted disks and unit-capacity points: we 
lift the problem to 3d, we construct an approximate conflict graph (as done by Mustafa 
and Ray), and we use a local search algorithm. This result also implies a PTAS for 
PackPoints for unweighted points and uniform capacity halfspaces in IR^. See Section [6] 
for the details. 

Hardness. We show some hardness results for our problems. In particular, we show that 
PackPoints for fat triangles in the plane is as hard as independent set in general graphs 
(see Lemma 17.21) . We also show that PackRegions is APX-hard (and thus there is 
no PTAS) for similarly sized fat triangles in the plane (thus "matching" the result of 
Corollary [ 



Main technical contribution. Besides the results mentioned above, our work further de- 
velops and extends the techniques for rounding LPs that rise out of low dimensional geometric 
constraints. Such work relies on finding the right order of making decisions about regions 
as they are being added, usually initially picking the elements to be considered randomly 
according to the value assigned to them by the associated LP. Such work in the context of 
LP rounding in geometric settings includes [VarlO| [CH091 ICCH09] . The basic idea is to build 
a conflict graph, on the appropriate random sample, and argue that there exists a vertex 
of low degree that can be added without throwing away too many conflicting vertices. Our 
work extends this approach to more involved settings where conflicts are not just whether 
two regions intersect or not (i.e., independent set in a graph), but rather involve a larger 
number of regions. To this end, we prove a combinatorial bound on the expected number of 
conflicts realized if we round the associated LP. A special easier case of this was addressed 
by Chan and Har-Peled |CH09j when the LP is an independent set LP. Naturally, in our case 
the analysis is considerably more involved. 

Previous work. Fox and Pach [FPU] presented an rf approximation for independent set 
for segments in the plane. The usage of LP relaxations for approximating such problems is 
becoming more popular. In particular, Chalermsook and Chuzhoy |CC09] use a natural LP 
relaxation to get an 0(log log m)-approximation for independent set of axis parallel rectangles 
in the plane. The geometric set cover problem and the more general problem, the geometric 
set multi-cover problem, have approximation algorithms that use e-nets to round the natural 
LP relaxation; see |CCH09j and references therein. Chan and Har-Peled |CH09] used local 
search to get a PTAS for independent set of pseudo-disks. Independently, Mustafa and Ray 
|MR10] used similar ideas to get a PTAS for hitting set of pseudo-disks in the plane. There 
is not much work on the hardness of optimization problems in the geometric settings we 
are interested in. |CC07] shows that the problem of independent set of axis-parallel boxes 
in three dimensions is APX-hard (the problem is known to be NP-Hard in the plane). See 
also |GCllt IHar09] and references therein for some recent hardness results. Naturally, in 
non-geometric settings, there is a vast literature on the problems and techniques we use, see 
|WSllj . Our algorithms use the randomized rounding with alteration technique to round 
a fractional solution. This technique was used in |Sri01j to find an approximate solution 
to packing integer programs (PIPs) of the form {maxwx : Ax < 6, x G Z'"}, where A is 
a matrix whose entries are either or 1. The approximation guarantee given in [SriOlj is 



0(n^/-^), where B = mirijfej. 

Organization. In Section [2] we define tlie problem and tlie associated LP relaxation, and 
describe some basic tools used throughout the paper. In Section [3] we present the approxi- 
mation algorithm for the hypergraph case. In Section H] we present various applications of 
our main result. In Section [5] we present the algorithm for packing points into fat triangles. 
In Section [6] we present a PTAS for some restricted cases. In Section [7] we present some 
hardness results. We conclude in Section |8l 

2 Preliminaries 

For a maximization problem, an algorithm provides an a -approximation if it outputs a 
solution of value at least opt/a, where opt is the value of the optimal solution. An (a,/3)- 
approximation algorithm for HGraphPacking is an algorithm that returns a (potentially 
infeasible) solution of value at least opt/a such that each hyperedge / contains at most 
max(7^(/), /3) vertices of the solution. 

a-fat triangles. For a > 1, a triangle A is a- fat if the ratio between its longest edge 
and its height on this edge is bounded by a (there are several equivalent definitions of this 
concept). A set of triangles is a- fat if all the triangles in the set are a-fat. The union 
complexity of n a- fat triangles is 0(nlog*n) [EASllt lAdBESllb] (the constant in the O 
depends on a, which is assumed to be a constant). 

2.1 LP Relaxation and the Rounding Scheme 

We consider the following natural LP relaxation for the HGraphPacking problem. For 
each vertex v, we have a variable x^ with the interpretation that x^ is 1 if f is selected, and 
otherwise. For each hyperedge /, we have a constraint that enforces that the number of 
vertices of / that are selected is at most the capacity of /. 

Hypergraph-LP : max 
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The energy of a subset X C V is £(X) = 2_, ^v- In the following, £ denotes the energy of 

the LP solution; that is £ = £(V) = Xlugv^^- Note that the energy is at most the number 
of vertices of the hypergraph. Also, we assume that £ > 1 (which is always true since all the 
capacities are at least one). 



Definition 2.1 Let G = (V, E) be a hypergraph. For any integer k, let -Ffc(') denote the 
function 

FAt)= max 1/ f e£ and\Xr]f\ = k + l\ 
xcy,\x\<t I J 

that is, Fk{t) is the maximum number of hyperedges of size k + 1 of a sub-hypergraph of G 
that is induced by a subset of at most t vertices. We say that G has bounded growth if 
the following conditions are satisfied: 

(A) There exists a non- decreasing function 7(- ) such that Fk{t) < 2'^^^H'y(t) for any k 
and t. 

(B) There exists a constant c such that Fk{xt) < cFk(t) for any t,k and x such that 
1 < a:< 2. 

This notion of bounded growth is a hereditary property of the hypergraph, and it is somewhat 
similar to the bounds on the size of set systems with bounded VC dimension. Hypergraphs 
with bounded growth arise naturally in geometric settings. 

The minimum capacity of a packing instance is a useful measure of how hard the instance 
is; formally, the minimum capacity of a given instance G is 

z/ = KG)=iinn#(/). (1) 

2.2 Basic tools 

Let G = (V, E) be a hypergraph, and let X C V be a subset of its vertices. The sub- 
hypergraph of G induced by X is G^ =(X, </nX /gEH. 

The following lemma testifies that the packing problem can be solved in a straightforward 
fashion if all the capacities are the same (i.e., uniform capacities). This is done by repeatedly 
applying a procedure to find and remove an independent set in the remaining induced sub- 
hypergraph (for example the procedure in jCHll] ). 

Lemma 2.2 Let G be a hypergraph for which there is a polynomial time algorithm alg that 
takes as input a fractional solution to Hypergraph-LP for an HGraphPacking instance 
on G, or any induced subgraph ofQ, with unit capacities — i.e., an independent set instance 
— and it constructs an integral solution whose value is at least an a fraction of the value 
of the fractional solution. Then one can compute in polynomial time a 2a -approximation 
for any instance o/ HGraphPacking on G with uniform capacities (i.e., all the hyperedges 
have the same capacity, say k). 

Proof: Let G = (V, E) be an instance of HGraphPacking in which all hyperedges have 
the same capacity, say k. Let Gq = G, and in the ith iteration, for i = 1, . . . ,k, compute 
a maximum weight independent set Yi in Gj„i using alg, and let Gj = Gx\Ui, where Ui = 
Yi U . . . U Fj. We claim that Uk is the required approximation. 

Clearly, no hyperedge of G contains more than k vertices of Uk as it is the union of k 
independent sets, and as such it is a valid solution. Now, let V"opt be the optimal solution. 
If w{Vopt n Uk) > w(yopt)/2 then we are done. Otherwise, consider the hypergraph Gi_i, 



and observe that V^pt \ f/j_i is a valid solution for HGraphPacking for this graph (with 
uniform capacities k). Interpreting this integral solution as a solution to the LP, and scaling 
it down by k, we get a factional solution to the independent set LP of this hypergraph of 
value w{Vopt \ Ui-i)/k. Since Yi was constructed using alg on the optimal fractional solution 
to the independent set LP of this hypergraph, we have that 

w{Yi) > optLp(Gi_i)/a > wiVopt \ Ui^i)/ka > wiVopt \ Uk)/ka > w(yopt)/'2ka. 

Which implies that w{Uk) > w{Vopt)/'2a. m 

A hypergraph G = (V, E) shatters X C V if the number of hyperedges in Gx is 2'"^'. 
The VC dimension of G is the size of the largest set of vertices it shatters. 

The following is a "sparsification" lemma. Here we get better bounds than the standard 
technique, as we are using stronger sampling results known for spaces with bounded VC 
dimension. 

Lemma 2.3 Let G = (V, E) be an instance o/ HGraphPacking with VC dimension d, and 
consider its fractional LP solution of value opt and with energy £. Then, one can compute, 
in polynomial time, a valid fractional solution for the LP of G such that: 

(A) The value of the new fractional solution is > opt/ 12. 

(B) The number of vertices with non-zero value is 0{dE log£). 

(C) The value of each non-zero variable is equal to i/M for some integer i < M, where 
M = 0((i log £). 

(D) The total energy in the new solution is 0(£). 

Proof -.Let e = l/£, where £ = "^^x^, and x^ is the value the LP assigns to f G V in the 
optimal LP solution. Let T = 0{d\ogE) and R be a random sample of V of (expected) size 
r = £T = 0{{d/e) log(l/5)), created by picking each vertex v independently with probability 
Xy ■ T\j. This sample is a relative (e, l/2)-approximation |Harlll fHSllj . with probability of 
failure < pi = e'^^'^\ That implies that for any hyperedge / G E such that x{f) = J2v&f ^y 
we have 

|Rn/l „ , ,,„,/x{f) 



|R| ^(i + V2)^ J 

To interpret the above, observe that E[|R fl /|] = x{f)T and E[|R|] = £T, as such, a rough 
estimate of the expectation of |Rn/|/|R| is x{f)/8,. Thus, the above states (somewhat 
opaquely) that no hyperedge is being over-sampled by R. 

Since the expected size of R is r, by Chernoff's inequality, we know that |R| < 2r with 
probability at least 1 — p2, where p2 = e'^^'^^) (as £ > 1). Now, consider a hyperedge / 
with capacity k, and observe that x{f) < k. As such, |R fl /| < (1 + l/2)(a;(/)/£ + e) |R| < 
(3/2) (/c + l)e2T < 6kT. In particular, if v appears t^ times in R (R is a multiset), then we 
assign it the fractional value y^ = t^/GT. We then have that y{f) = J2vef Vv — \^^ f\ /6^ ^ 



^A minor technicality is that XyT might be larger than one. In this case, we put [xijTj copies of v into 
R, and we put an extra copy of v into R with probability XyT — \_XyT\ . It is straightforward to verify that 
our argumentation goes through in this case. Observe that such large values work in our favor by decreasing 
the probability of failure. 



k (and this holds for all hyperedges with probability > 1— pi). As such, the fractional solution 
defined by the |/'s is valid. 

As for the value of this fractional solution, consider the random variable Z = ^^ y^w^v), 
which is a function of the random sample R. Observe that 
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as optLp = J2v'^i'^')^v- In particular, since no vertex can have w{v) > optLP (otherwise, we 
would choose it as the solution), it follows that 

V — ^ / X V~~^ ty R 2t 2£T C 

Z = 2^yvW{v) < optLp 2^-^< OptLP -^ < OPtLPTT^ = ^pt^p-r;^ < -OptLp- 



This implies that Pv[Z > optLp/12] > l/4£ = e/4. Indeed, if not. 
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a contradiction. As such, a random sample R corresponds to a valid assignment with value 
at least optLp/12 with probability at least Pr[Z > optLp/12] — pi — P2 > s/8, as pi + p2 
is an upper bound on the sample R failing to have the desired properties. As such, taking 
u = 0(£log£) independent random samples one of them is the required assignment, with 
probability > 1 — (1 — e/8)" > 1 — l/£*^''^^ We take this good sample together with its 
associated LP values as the desired fractional solution to the LP. Also, note that the total 
energy of the new solution is 6(£), since by Chernoff's inequality r/2 < |R| < 2r with 
probability at least 1 — 2p2. ■ 

3 Approximate packing for hypergraphs 

In this section, we present the algorithm for computing a packing for a given hypergraph 
G = (V, E). We assume that |E| is polynomial in |V| and that G has the bounded growth 
property introduced in Definition 12.11 (properties which both hold for hypergraphs arising 
out natural geometric settings). Let x be a solution to the Hypergraph-LP relaxation 
described in Section [2711 

3.1 The algorithm 

We round the fractional solution to an integral solution using a standard randomized round- 
ing with alteration approach. The first step is to choose an appropriate ordering of the 
vertices. We will see later how to choose a good ordering; for now, we assume that we are 
given the ordering. The rounding then proceeds in two phases, the selection phase and 
the alteration phase. In the selection phase, we pick a random sample C of the vertices 
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by picking each vertex v independently at random with probabihty Xj,/A, where A is a 
parameter that we will determine later. In the alteration phase, we pick a subset of C as 
follows: We consider the sampled vertices in the order chosen and we add the current vertex 
to our solution if the resulting solution remains feasible. We say that a vertex is selected 
if it is present in the sample, and we say that it is accepted if it is present in the solution. 
The main insight is that we can take advantage of the bounded growth property of the hy- 
pergraph to show that there is an ordering such that each vertex is accepted with constant 
probability, provided that it is selected. This will immediately imply that the algorithm 
achieves a 0(A)-approximation. 

The main challenge is to prove that a good ordering for the alteration phase exists, that 
is an ordering such that we accept each selected vertex with constant probability. We now 
proceed to give such a proof. This proof will suggest a natural 0{'nP^^^^') time brute force 
algorithm to actually compute this good ordering, where C is the maximum capacity of an 
edge in the given instance. In Section 13.31 we show how one can remove this exponential 
dependence on C. 

Running Example 3.1 To keep the presentation accessible, we interpret this algorithm 
for instances of PackRegions in which the regions are disks. Specifically, we are given a 
weighted set of disks 2) and set of points P with capacities. The disks of T> form the set of 
vertices of the hypergraph and every point p G P forms a hyperedge; that is the hyperedge /p 
is the set of all disks of 1) that contain p. 

In this case, the mysterious quantity Fk{t) (see Definition \2.1\) is hounded by the number 
of faces in an arrangement of t disks that are of depth exactly k + 1. Since the union 
complexity of t disks is linear, standard application of the Clarkson technique implies that 
Fk{t) = 0{kt). In particular, for this case 7(t) = 0(1). 

3.1.1 Constructing a good ordering 

Before we describe how to construct a good ordering of the vertices, it is useful to understand 
what will force a vertex to be rejected in the alteration phase. With this goal in mind, 
consider an ordering of the vertices. Let C be a sample of the vertices in V such that each 
vertex w is in C independently at random with probability x^/A. Let v be a vertex in C. 
When we consider v in the alteration phase, we will reject v iff there exists a hyperedge / 
of capacity #(/) such that / contains v and we have already accepted #(/) vertices of /. 
The event that we already accepted #(/) vertices of / is difficult to analyze. However, as 
we will see, we can settle for a more conservative analysis that upper bounds the probability 
that V is rejected, given that all of the vertices in C that appear before v in the ordering are 
accepted. (In the alteration phase, it is possible that not all vertices in C that appear before 
V will be accepted, but this can only help us.) Since we are only interested in the event that 
C contains k + 1 vertices — the vertex v and k other vertices that appear before v in the 
ordering — that are contained in a hyperedge of capacity k, only the set of vertices that 
appear before v in the sample matter, and not the actual ordering of the vertices. With this 
observation in mind, we define a k-conflict to be a set of /c + 1 vertices that are contained 
in a hyperedge of capacity k. In the following, Hk denotes the set of all /c-confiicts, and 
Ti = UkHk denotes the set of all conflicts. We are interested in the probability of the event 
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that all of the vertices of a /c-conflict, /i, are present in the sample, and we refer to this 
probability as the A-potential of the conflict, p/^{h). For the analysis it will also be useful 
to define the unsealed version of this quantity, that is the probability that all the vertices of 
a conflict are present given that we sampled each vertex with probability x^ and not Xt,/A. 
We refer to this quantity as simply the potential of the conflict, p{h). Formally we have. 



PA{h) = n ^ ^^^ p^^^ = n 



Another quantity of interest is the expected number of conflicts in which a vertex v partic- 
ipates, given that v is in the sample. We refer to this quantity as the A-resistance of a 
vertex w in a set of vertices X C V, and we use r]A{v, X) to denote it: 

r]A{v,X) = — V pA{h). 

'" h&H,hCX,veh 

The ordering. Note that, if the A-resistance of v with respect to the set X of vertices 
that come before it in the ordering is small, the probability of rejecting v is also small. This 
suggests that the vertex with least resistance (with respect to V) should be the last vertex 
in the ordering. This gives us the following algorithm for constructing an ordering: We 
compute the vertex of least resistance and put it last in our ordering (i.e., it is Vn)- We then 
recursively consider the remaining vertices and we compute an ordering for them. In the 
following, we assume for simplicity that the resulting ordering is Vi, . . . ,u„. 

Note that computing the resistance of a vertex by brute force takes ©(n'"^*^'-^^) time, 
where C is the maximum capacity of a hyperedge, and therefore this algorithm is not efficient. 
We give a polynomial time algorithm for constructing the ordering in Section 13.31 

3.2 Analysis 

Our main insight is that, if the hypergraph satisfies the bounded growth property defined in 
Definition 12. H then for any set X C V there exists a vertex f G X such that ?7a(^, X) < 1/4. 
We prove this below in Section 13.2.21 (see Lemma 13. 9p . This proof requires that we set 
A = a7(£)^/'^, where a is some sufficiently large constant. As such, in the remainder of this 
section we assume A = a7(£)^/'^. 

We now show that given ?7a(w,X) < 1/4, proving the quality of approximation of the 
algorithm is straightforward. 

Lemma 3.2 Let C and O he the set of vertices that were selected and accepted by the algo- 



rithm, respectively. For each i, we have Pr 



ViGO 



Vi G C 



>3/4. 



Proof: Let Xi = {vi, . . . ,Vi). Note that, if we selected Vi, we rejected Vi in the alteration 
phase only if Vi participates in a confiict with some of the vertices in {vi, . . . , f j_i} fl C. Let 
Zi be the number of confiicts of Xi that contain Vi and are realized in C, i.e., h 'O C. In the 
following, we show that the probability that Z^ is non-zero is at most 1/4, which implies the 
lemma. 
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Consider a /c-conflict h = {vj^, . . . , Vj^,Vi}, where each vertex of h is in Xi and h contains 
Vi. The probabihty that all of the vertices of h are selected, given that Vi is selected, is equal 



to -^ ■ -^ ^ = — PAih). Therefore we have 

A A A ~ 



E 



X7' 



Zi 



ViEC 



E 



heH,hCXi,Vieh 



A 1 

—pA{h) = VA(,Vi,Xi) < -, 

x, 4 



where the last inequality follows from Lemma 13.91 and the fact that Vi is the vertex of 
minimum resistance in X^. Thus 



Pr 



V^^O 



Vi e C 



<Pr 



Zi > 



Vi e C 



<E 



Vi e C 



< 



Therefore, if Vi is selected, it is accepted with probabihty at least 3/4. ■ 

Corollary 3.3 The total expected weight of the set of vertices output by the algorithm is 
r2(opt/7(£)^/'^), where opt is the weight of the optimal solution, and v is the minimum 
capacity of the given instance. 



Proof: By Lemma [3.21 for each vertex f G V, we have 



Pr[-y eO]= Pr[(t; e O) n{v e C)] = Pr 



veO 



V eC 



Pr[t; eC]> -Pr[t; G C] 



> 



4A 



where A = 0(7(£)^/''). By linearity of expectation, we have that the expected weight of the 
generated solution is at least 



dGV 



4A 



Wo, = Q 



7 y ll X^W^, 

7(£)i/- 



Q 



opt 



7(6)1/-; ' 



as '^^XvWv is the value of the fractional LP solution, which is bigger than (or equal to) the 
weight of the optimal solution. ■ 

3.2.1 On the expected number of conflicts being realized 

To analyze the algorithm we need to understand how conflicts might form during its exe- 
cution, and show that the damage of such conflicts to the generated solution is limited. To 
this end, consider the quantity 



Ffc(t) = max 

-4CX,|A|<i 



{/|/GEand \Anf\ = k + l] 



This is the maximum number of /c-confiicts that can be realized for a set of t vertices. The 
quantity of interest in the following is 'YIiK^'h Pi^)^ ^ ^^ ^^ ^^^ expected number of conflicts 
that would be realized if we sample according to the LP solution. Our purpose is to prove 
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that this quantity is bounded by a function of the energy of the LP (the bound would involve 
the function -Ffc(-) defined above). 

With this goal in mind, we let R be a random sample of X such that each vertex v & X 
is in R independently at random with probability Xy/2. We stress that R is a random sample 
that we use for the purposes of defining a quantity M (i.e., the expected number of conflicts 
realized in R), and it should not be confused with the random sample C that is used by 
the algorithm. In the following, we bound M from above in Lemma 13.41 and from below in 
Lemma [3.51 Putting these two bounds together imply the desired bound on Ylheu Pi^)- 

A conflict /i G "H is realized in R if there is a hyperedge / € E such that h = f r\R and 
\h\ = ifif) + l. 

The following is similar in spirit to the Clarkson technique (a similar but simpler argument 
was used by Chan and Har-Peled |CHllj ). 

Lemma 3.4 The expected number of k- conflicts realized in R is M. = 0(-Ffc(£(X))), where 
R is a random sample of X such that each vertex v & X is in R independently at random 
with probability x^/2. 

Proof: Each /c-conflict h that is realized corresponds to a hyperedge / with capacity k such 
that h = f nR. Additionally, two realized conflicts that are distinct correspond to different 
hyperedges. Therefore the number of A;-conflicts that are realized in R is at most the number 
of hyperedges / such that the capacity oi f is k and \f r\R\ = k + 1. It follows from the 
definition of -Ffc(-) that the number of /^-conflicts is at most -Ffc(|R|). Therefore it suffices to 
upper bound the expected value of Ffc(|R|). 



Note that E[|R|] = £(X)/2. We have 



E 



F,(|R|)] <£ Pr[|R|>t^]F.((t + l)^) <f;2-(*+^)/^F,((t + l)^) 



oo 

^ ^2-(*+i)/2cOa°g*)F4£(X)) = 0(Ffc(£(X)) 



t=o 

since G has the bounded growth property (see Definition 12. ip . and by the Chernoff inequality 
(we use here implicitly that £(X) > 1). ■ 

Lemma 3.5 For each k-conflict h, the probability that h is realized in R is at least p{h) /2(2e)'^ 
Therefore the expected number ofk-conflicts realized inRisJA = Vtii '^h&n hex P(^) ) /(^s)'^ 

Proof: Let / G E be a hyperedge with capacity k that generated the confiict h. Since x is 
a feasible solution for the LP, we have that ^^^f_,^Xy < J2vef ^^ — i^if) — ^- Clearly, the 
confiict h is realized if we pick all the vertices of h, and none of the vertices oi f — h, and 
the probability of that event is 

p{h) 



2 J-i V 2 / - 2^+ 

vGh v£f—h vGh 



■ exp 




2fc+i -^^^\ ^ -- - 2(2e)fc' 
In the first line we used the inequality 1 — a;„/2 > exp(— a;„), which holds since Xy < 1. 
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Putting the above two lemmas together, we get the following. 

Lemma 3.6 For any non-negative integer k we have N, P(^) — 0[{2e) Ffc(£(X)) 

heUkAQX 

Running Example 3.7 In our running example, we have that the expected number of k- 
conflicts that are being realized by a random sample (sampling more or less according to the 

LP values) is Ylhen P(^) ~ 0((2e)^/c£]. This is a hefty quantity, but the key observation 
is that if we sample according to the LP values scaled down by a large enough constant, then 
the probability of such a conflict to be realized drops exponentially with k. In particular, for a 
sufficiently large constant, the expected number of realized k-conflicts in such a smple is going 
to be < £/ (lO ■ 2^^) . Intuitively, this implies that such conflicts can only cause the algorithm 
to drop very few vertices during the rounding stage, thus guaranteeing a good solution. 

3.2.2 Resistance is futile, if you pick the right vertex 

In the following, we consider a subset X of the vertices and we show that there exists a 
vertex v E X whose A-resistance ri^{v,X) is at most 1/4. Recall that Tik is the set of all 
fc-conflicts involving vertices in V. We can rewrite the A-resistance of f in X as 

A 1 °° 1 

^ hG'H,hCX,vGh " k=u h<^Hk,hCX,v£h 



As shown in Lemma 13.61 we can relate the total potential of the conflicts of Tik that are 
contained in X to the maximum number of /c-confiicts contained in a set of at most 8,{X) 
vertices, where £(-^) = J2vex-'^^- 

Recall that the hypergraph has the bounded growth property (see Definition 12. ip and 
this property is hereditary. Therefore the function -Ffc(-) in the lemma above has the two 
properties described in Definition 12.11 and we get the following corollary. 

Corollary 3.8 We have Y^hen hex Pi^) ^ 0(2'^'''£(X)7(£(X))), where c is a constant. 
We can use Corollary 13.81 to complete the proof of Lemma 13.91 as follows. 

Lemma 3.9 Suppose that the hypergraph G satisfies the bounded growth property (see Def- 
inition lKl\) . Let A = a'y{8,y^'^, where a > is a sufficiently large constant and v is the 
minimum capacity of the given instance (see Eq. (1\)). Then, for any set X C V, there exists 
a vertex v E X such that ri^{v,X) < 1/4. 

Proof: Let T = J2vgx •'^'"Va{v,X). The quantity T/£(X) is the weighted average of the 
resistances of the vertices in X, where the weight of a vertex v is Xv/8,{X). Therefore it 
suffices to show that T < £(X)/4, since the minimum resistance is at most the weighted 
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average. We have 

k=u vex heUk k=y heUk k=u h&Hk 

hex hex hex 

v£h 

k=u k=u 

"■ V ' 

=5 

by Corollary 13 -St where (3 is some constant. Since A = a7(£)^/'^, we have 

°° /OCX fc / 1 \ '^/'^ °° /r)c\k / ^ \k/u -. 

^-eK^)(^-<^) ^(^w)^E/'y (^-<^) ^(e)^j. 

In the second to last inequality, we have used the fact that 7(- ) is non- decreasing. The last 
inequality follows if we pick a to be a sufficiently large constant. Therefore T < £(X)/4, 
and the lemma follows. ■ 

3.3 Improving the running time 

In Section |3l we described an algorithm that constructs an ordering of the vertices by re- 
peatedly finding the vertex of least resistance with respect to the set of remaining vertices. 
Computing the resistance of a vertex by brute force takes 0{n'~^~^'^^^^) time, where C is the 
maximum capacity of an edge in E. However, for our analysis to go through, we only need 
to find a vertex that is safe with respect to the set of remaining vertices; informally, a vertex 
V is safe if the probability that it participates in a confiict with a random sample of the 
remaining vertices is smaller than some constant (that is strictly smaller than one), where 
each remaining vertex u is included in the sample with probability x^j/A. In this section we 
show that there is a sampling algorithm that finds a safe vertex with high probability and 
its running time is polynomial in the maximum capacity C. 

Lemma 3.10 Computing a good ordering of the vertices can be done in polynomial time. 
Namely, the algorithm of Section\^can he implemented in polynomial time. 

Proof: To get the same quality of approximation we do not need to take the vertex of least 
resistance in each round (of computing the ordering), but merely a vertex that is "safe." 
More precisely, let X be the current set of vertices, let w be a vertex of this set, and let R 
be a random sample of X in which each vertex u is included with probability Xu/ ^ (also we 
force V to be in R). We say that v is violated in R if f is contained in a hyperedge / such 
that the number of vertices of / that are in R is larger than its capacity #(/). Let fi{v,X) 
denote the probability that v is violated in R. Note that fi{v,X) is a (conservative) upper 
bound on the probability that v is rejected by our rounding algorithm if we started with an 
ordering in which X \ {v} is the set of all vertices that come before v. Therefore, in order 
for our rounding to succeed, in each round we only need to find a vertex v for which the 
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probability fi{v,X) is low, where X is the set of all vertices that still need to be ordered at 
the beginning of the round. (We remark that it follows from the argument of Lemma 13.91 
that, for any set X, there is a vertex v for which fi{v,X) < 1/4.) 

Now we are ready to describe how to construct an ordering for our algorithm. Let X 
be the set of vertices that still need to be ordered. As we will see shortly, for each vertex 

V ^ X, we can compute an estimate /l(w,X) of the probability fi{v,X). We pick the vertex 

V with minimum estimated probability fl^v.X), we make v the last vertex (in the ordering 
of X) and we recursively order X \ {u}. 

We can compute the estimates 'Jl{v^X) in polynomial time as follows. Fix a vertex v. 
Let iphe a, sufficiently large polynomial in n. We pick ip independent random samples of X 
(again, forcing v to be in each of these samples); in each random sample, each vertex u is 
included with probability x„/A. We set Jliv.X) to be the fraction of the samples in which 
the vertex v is violated. Using a standard argument based on the Chernoff inequality, we 
can show that our estimates are very close with high probability, and therefore our rounding 
algorithm achieves the required approximation with high probability as well; we omit the 
easy but tedious details. ■ 

3.4 The result 

Theorem 3.11 Let G = (V, E) he a hypergraph with a weight function w{-) on the vertices 
and a capacity function 7^(-) on the edges, such that |E| is polynomial in |V| and G satisfies 
the bounded growth property (see Definition \KT\); that is, Fk{t) = 2'^''^'t'y(t). Then we can 
compute in polynomial time a subset X O \/ of vertices such that no hyperedge f contains 
more than its capacity i^{f) vertices of X. Furthermore, in expectation, the total weight of 
the output set is n(opt/7(£)^/''), where opt is the weight of the optimal solution, and v is 
the minimum capacity of the given instance. 

Consider an integer constant > 0, and observe that one can always relax the capacity 
constraints of a given instance of HGraphPacking by replacing all capacities smaller than 
(f) by 0. Theorem 13 . 1 1 1 thus implies the following. 

Corollary 3.12 Given an instance of HGraphPacking, with the bounded growth prop- 
erty, one can compute in polynomial time a (0(7(6)^/'^) , (fj -approximation to the optimal 
solution. 

3.5 Contention resolution schemes 

Chekuri et al. jCVZll] considered a broad class of rounding schemes, which they called 
contention resolution schemes (CR schemes). Informally, the family of all CR schemes 
consists of all rounding strategies based on randomized rounding with alteration. The precise 
definition of a CR scheme is the following. 

Let A^ be a finite ground set of size n, and let / : 2^ — t- IR+ be a submodulacl set 
function over A^. Let I C 2^ be a downward-closeq^ family of subsets of A^. The problem of 



^A function / : 2^ -)• M is submodular if f{A) + f{B) > f{A n B) + f{A U B) for any two subsets A, B 
of N. Additionally, / is monotone if f{A) < f{B) for all subsets A, B such that A C B. 
^A family I of subsets of TV is downward-closed if i? e I and A (- B then A £l. 
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maximizing f{X) subject to the constraint that X G X generahzes the hypergraph packing 
problem: the function / satisfies f{X) = J2vex ^v, and the family X is the family of all 
subsets X C y such that, for each hyperedge /, |/ fl X| < #(/). Let Pi C [0, 1]" be a 
convex relaxatioio of the constraints imposed by X; the set of all feasible fractional solutions 
to the Hypergraph-LP relaxation described in Section [2TT] is a convex relaxation for the 
family of all feasible solutions to the hypergraph packing problem. Let F be the multilinear 
extensiorl^ of /. Let x be a feasible solution to the relaxation {maxF(x) : x G Px}- The 
definition of the multilinear extension F suggests the following natural rounding strategy: 
given X, we construct a random set R{x) by picking each i & N independently at random 
with probability Xj. The expected value of f{R{x)) is equal to F{x), but it is unlikely that 
R{x) is in X. To address this, we want to remove some elements from R{x) in order to get 
a subset / C R[x) such that I EX. We want this step to have the property that, for each 
i E N, the probability that i is in / is at least cXi, for some parameter c > 0. Chekuri et al. 
[CVZllj call such a rounding strategy a c-balanced CR scheme for Px- In certain settings it 
is convenient to scale the fractional solution; the rounding strategy described above for the 
hypergraph packing problem is one such example. This motivates the following more general 
CR scheme. 

Definition 3.13 ( |CVZlT| ) A {b,c) -balanced CR scheme for Px is a scheme such that for 
any x G Px, the scheme selects an independent subset I C R{bx) with the following property: 
Pr[i G / I i G R{bx)] > c for every element i E N. The scheme is said to be monotone if 
Pr[z G / I R{bx) = Ri] > Pr[z G / | R{bx) = R2] whenever z G -Ri C i?2- A scheme is said to 
be strict z/Pr[z G / | z G R{bx)] = c for every i. 

Chekuri et al. showed that, if / is the output of a monotone [b, c)-balanced CR scheme, the 
expected value of / is at least cElF^bx)]. 

Theorem 3.14 Let / : 2^ — )■ IR+ be a non-negative submodular function and let x be a 
point in Px, where Px is a convex relaxation fori C 2^ . Let I{x) eX be the random output 
of a monotone {b,c) -balanced CR scheme on x E Px- If f is non-monotone, let us assume 
in addition that the CR scheme is strict. Then E[/(/(x))] > c'E[F{bx)]. 

The rounding scheme described in Section [3?T] is a monotone (A, l/4)-balanced CR scheme 
on X G Px, where Px is the set of all feasible solutions to the Hypergraph-LP relaxation. 
Therefore Theorem 13.111 extends to the setting in which the total weight of a set of vertices 
is a monotone submodular function instead of a linear function. 

Corollary 3.15 Let G = (V, E) be a hypergraph. Let w : 2^ ^ IR+ be a weight function on 
the vertices that is non-negative, monotone, and submodular. Let #(•) be a capacity function 
on the edges such that |E| is polynomial in |V| and G satisfies the bounded growth property (see 
Definition \2.1\) : that is, Fk{t) = 2^^''H'j{t). Then we can compute in polynomial time a subset 
X (^ y of vertices such that no hyperedge f contains more than its capacity #(/) vertices 
of X . Furthermore, in expectation, the total weight of the output set is f2(opt/7(£)^''^), 



^°Pi is the closure of the set of characteristic vectors of the sets in I under convex combinations. 

"The muhihnear extension F : [0,1] 1^1 ^ M of a function / : 2^ ^ R is the function F{x) = 
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where opt is the weight of the optimal solution, and v is the minimum capacity of the given 
instance. 

4 Applications 

Using our main result (Theorem I3.1ip . we get several approximation algorithms for the 
packing problems mentioned in the introduction. We present some of these results here. 

4.1 Packing regions with low union complexity 

Let 2) be a set of n weighted regions in the plane, and let the maximum union complexity 
oi m < n objects of T> be U(?7i) = ?7iu(m). We assume that (i) U(n) /n = u(n) is a non- 
decreasing function, and (ii) there exists a constant c, such that U(xr) < cU(r), for any r 
and 1 < X < 2. We are also given a set of points P, where each point p G P is assigned a 
positive integer #(p) which is the capacity of p. 

We are interested in solving PackRegions (Problem II. ip for T) and P. Consider the 
hypergraph G obtained by creating a vertex for each region and a hyperedge for each subset 
of regions containing a given point of P. Here, Fk{t) is bounded by the number of faces in the 
arrangement of t regions of depth exactly k + 1. The number of such faces can be bounded 
by the union complexity by a standard application of the Clarkson technique [Cla88[[CS89j . 



Lemma 4.1 Consider a set of regions in the plane such that the boundary of every pair 
intersects a constant number of times. The number of faces of depth at most k + 1 in the 
arrangement of any subset of these regions of size t is 0{k'^\J{t/k)). 

Plugging this bound into Theorem 13.111 yields the following result. 

Theorem 4.2 Let T) be a set of m weighted regions in the plane such that the union com- 
plexity of any t of them is U(t) = tu{t). Let P be a set of n points in the plane, where there 
is a capacity #(p) associated with each point p G P. There is a polynomial time algorithm 
that computes a subset O ^ T) of regions such that no point p E P is contained in more 
than #(p) regions of O. Furthermore, in expectation, the total weight of the output set is 

r2(opt/u(£) ], where opt is the weight of the optimal solution, £ is the energy of the LP 
solution and v is the minimum capacity of the given instance. 

Alternatively, for any integer constant (j), one can get a iO { u(£) ' j , j -approximation 
to the optimal solution for the given instance. 

The following results follow from the theorem above. 

Corollary 4.3 (A) The union complexity of pseudo-disks and fat triangles of similar size is 
linear; that is, U(t) = 0{t). Therefore we get an 0{1)- approximation for PackRegions if 
the regions are fat triangles of similar size, disks, or pseudo-disks. 

(B) The union complexity of fat triangles is 0(n log* n) lAdBESllaf . Therefore we get 
an O (log* n) - approximation for PackRegions if the regions are (arbitrary) fat triangles. 

(C) Consider a set of regions in the plane such that any pair of them intersects a constant 
number of times (e.g., a set of arbitrary triangles). In this case, U(t) = O(t^) and u(t) = 
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0{t). Therefore, for any integer constant (p > 0, we get an (©(S^/*^) , (fj -approximation for 
instances of PackRegions on such regions. 

4.2 Packing halfspaces, rays and disks 

Problem 4.4 (A) PackHalfspaces; Given a weighted set of halfspaces S and a set of 
points P with capacities in H^, find a maximum weight subset O of § so that, for each point 
p, the number of halfspaces of O that contains p is at most #(p). 

(B) PackRaysInPlanes; Given a weighted set of vertical rays "Ji and a set of planes "K 
with capacities in ]R^, find a maximum weight subset O of "Jl so that, for each plane h, the 
number of rays of O that intersect h is at most #(h). 

(G) PackPointsInDisks; Given a set T) of disks with capacities and a weighted set P 
of points, find a maximum weight subset O of the points so that each disk r & T) contains at 
most #(r) points of O. 

Since the union complexity of halfspaces in three dimensions is linear, we get the following 
from Theorem 13.111 (and the 3d analogue of Lemma [4. ip . 

Corollary 4.5 One can compute, in polynomial time, a constant factor approximation to 
the optimal solution of the PackHalfspaces problem. 

Standard point/plane duality implies that the same result holds for the dual problem. 
Namely, a point (a, 6, c) gets mapped to the plane z = ax + hy — c and a plane z = ax + by + c 
gets mapped to the point {a,b,—c). Also, a point lies below a given plane if and only if 
the dual point of the plane lies below the dual plane of the point. As such, the dual of an 
instance of PackHalfspaces is an instance of PackRaysInPlanes (and vice versa). 

Thus, Corollary 14.51 implies the following. 

Corollary 4.6 One can compute, in polynomial time, a constant factor approximation to 
the optimal solution of the PackRaysInPlanes problem. 

Finally, observe that an instance of PackPointsInDisks can be lifted into an instance 
of PackRaysInPlanes, by the standard lifting f{x, y) = {x, y, x"^ + y'^), which maps points 
and disks in the plane to halfspaces and points in three dimensions [dBCvKOOSj . 

Corollary 4.7 One can compute, in polynomial time, a constant factor approximation to 
the optimal solution to the PackPointsInDisks problem. 

4.3 Axis Parallel Rectangles/Boxes 

4.3.1 Packing rectangles (2d) 

Problem 4.8 (PackRectsInPoints.) Given a weighted set S of axis-parallel rectangles 
in the plane, and a point set P with capacities, find a maximum weight subset O (^ "B, such 
that, for any p G P, the number of rectangles of O containing p is at most #(p). 
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Note that the union complexity of a set of rectangles can be quadratic. Hence we cannot 
simply use Theorem 14. 21 to get a meaningful approximation. However, by using the standard 
approach for approximating the independent set of rectangles, one can get a reasonable 
approximation as the following lemma testifies. 

Lemma 4.9 Given an instance (S, P) of PackRectsInPoints with m rectangles, one 
can compute, in polynomial time, a subset O ^ "B of total weight r2(opt/log?Ti) such that no 
capacity constraint of P is violated, where opt is the weight of the optimal solution. 

Proof: It is straightforward to verify that a set of rectangles that intersect a common line 
have linear union complexity. Therefore it follows from Theorem 14.21 that we can get a 
constant factor approximation for such instances. Given an arbitrary set of axis-parallel 
rectangles, we can reduce it to the case in which all rectangles intersect a common line as 
follows. 

We construct an interval tree on !B. Let b e S be the median rectangle of S when 
sorted by left edges. Let £ denote the vertical line which passes through the left edge of b, 
and let "Bi denote the set of rectangles it intersects. We associate "Bi with the root of our 
tree, and then recursively build left and right subtrees for the rectangles in 23 \ 23^ that lie 
to the left or right of i, respectively. The recursion bottoms out once every rectangle has 
been stabbed by a line. Clearly the depth of the tree is 0(log?n), since each time we choose 
the median line and only recursively continue on those rectangles that do not intersect it. 
Therefore there exists a level of the range tree that has a solution of weight r2(opt/log?7i). 
The algorithm now considers each level of the tree separately, and for each node at the 
given level, it constructs an approximate solution for the rectangles associated with the node 
using the constant factor approximation algorithm guaranteed by Theorem 14. 2[ Next, the 
algorithm considers the union of all these solutions to form the solution for this level. Since 
two rectangles associated with two different nodes at the same depth in the range tree do 
not intersect, the resulting set is a valid solution. 

The algorithm returns the best solution found among all the levels. Clearly, its weight is 
0(opt/logm). ■ 

4.3.2 Packing axis-parallel boxes (3d) 

The union complexity of axis-parallel boxes in IR that contain a common point is also linear, 
and therefore a similar approach as above will enable us to solve the following problem. 

Problem 4.10 (PackBoxesInPoints.) Given a weighted set "B of axis-parallel boxes in 
]R^, and a point set P with capacities, find a maximum weight subset O C B, such that, for 
any p G P, the number of boxes of O containing p is at most #(p). 

Lemma 4.11 Given an instance (S, P) of PackBoxesInPoints with m boxes, one can 
compute, in polynomial time, a subset O '^ B of total weight r2(opt/log m) such that no 
capacity constraint of P is violated, where opt is the weight of the optimal solution. 

Proof: We build a multi-layer interval tree on the boxes. On the top layer, we build a 
balanced tree on the x-axis projection of the boxes, where a node Vx' stores all boxes inter- 
secting the plane x = x'. Next, we build for each such node a secondary interval tree on the 
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y-axis projections, and for each node on this secondary data-structure we build a third layer 
data-structure on the 2;-axis projections. All the boxes are stored in the nodes of the third 
layer data-structure. 

First, observe that the boxes stored in a node on the third layer, Vx'y'z'i all contain the 
point (x', y', z'). Now, the union complexity of axis parallel boxes all sharing a common point 
is linear. As such, we can apply Theorem 13.111 to compute a packing that does not violate 
the capacities and is a constant factor approximation to the optimal (on this restricted set 
of boxes). Now, for each third layer tree, we find the level that contains the best possible 
combined solution (taking the union of the solutions at a level is valid since there is no point 
in common to two boxes that are stored in two different nodes at the same level). This 
assigns values for each node in the secondary tree. Again, we choose for each secondary tree 
the level with the maximum total weight solution. We assign this value to the corresponding 
node in the first layer data-structure. Again, we choose the level with the highest possible 
value. This corresponds to a valid solution that complies with the capacity constraints. 

As for the quality of approximation, observe that for a tree at a given layer, at least a 
logarithmic factor of the remaining weight of the optimal solution is contained in some level 
of the tree. As such, each time we go down one layer in this data-structure we lose at most 
a logarithmic factor of the optimal solution, and hence the quality of approximation of this 
algorithm is fi(log m). ■ 

Remark 4.12 (A) It is natural to ask if this result can be extended to higher dimensions. 
However, it is easy to see that the union complexity of n axis-parallel boxes in four dimen- 
sions that all contain a common point can be quadratic. As such, this approach would fail 
miserably. We leave the problem of getting a better approximation for this case as an open 
problem for further research. 

(B) The continuous version seems to be considerably easier. For the weighted case (with 
unit capacities; that is, the independent set variant) an 0(log "^ m/ log log m) approximation 
is known [CH09J by solving the two dimensional case, and then using the above interval 
tree technique to apply it for higher dimensions. For the unweighted continuous case, an 
O (log log m log ~^m) approximation is known |CC09] . 

The discrete version is different than the continuous version because, for example, if 
considering boxes in IR^ that all intersect the x|/-plane, the induced two dimensional instance 
fails to encode the capacity constraints, as they rises from points in three dimensions that do 
not lie on this plane (while in the continuous case, it is enough to solve the induced problem 
in this plane). 

4.3.3 Packing points into rectangles 

Problem 4.13 (PackPntsInRects.) Given a weighted set P of points and a set S of 
axis-parallel rectangles with capacities in the plane, find a maximum weight subset O C P, 
such that, for any b G S, the number of points of O contained in b is at most #(b). 

We first observe that the hypergraph that arises from a instance G = (P, 3) of PackP- 
ntsInRects, might not have the bounded growth property for any reasonable growth func- 
tion. To see this consider the following example. 
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Consider two parallel lines in the plane with posi- 
tive slope. Place n/2 points on each line such that all 
the points on the top line lie above and to the left of 
all the points on the bottom line. Let the set of rect- 
angles for this instance of PackPntsInRects be all 
the rectangles which have a point on the top line as 
their upper left corner and a point on the bottom line 
as their lower right corner. In this case any subset of 
0{t) points from the top line and 0{t) points from 
the bottom line induce a set of O(t^) hyperedges, each 
of size 2. Therefore, Fi{t) = fi(t^), and hence Theo- 
rem [3lTT] only gives an 0(£) approximation. 

Since we cannot hope to apply our main result to the case of PackPntsInRects, we will 
instead seek a bi-criterion approximation. Our algorithm here is inspired by the work of Ezra 
et al. JAESIO] on e-nets for rectangles. Before tackling this problem, we will first consider 
an easier variant, which will be useful later in obtaining a bi-criterion approximation. In the 
following, we call a set of rectangles such that all their (say) bottom edges lies on a common 
line a skyline. 

Problem 4.14 (PackPntsInSkyline.) Given a weighted set P of points and a set S of 
skyline rectangles with capacities in the plane, find a maximum weight subset O C P, such 
that, for any b G S, the number of points of O contained in h is at most #(b). 

Lemma 4.15 Let P be a set of n points in the plane all placed above the x-axis. Let Fk{n) 
be the maximum number of different subsets of P of size k that are realized by intersecting P 
with a rectangle whose bottom edge lies on the x-axis. We have that Fk{n) = 0{nk'^). 

Proof: Consider a rectangle b with its bottom edge lying on the x-axis, and which contains 
k points of P. Lower its top edge till it passes through a point of P, and let p denote this 
point. Similarly, move its left and right edges till they pass through points of P. Let b' be 
this new canonical rectangle. Now, let ii^ft (resp. Zright) be the number of points of P inside b' 
that are to the left (resp. right) of p. Clearly, (p, zioft, bright) uniquely identifies this canonical 
rectangle. This implies the claim as p G P, iictt < k and Zj-ight < k, and hence the numbers of 
such triples is 0(|P| k^). m 

Lemma 4.16 Given an instance o/PackPntsInSkyline, one can compute, in polynomial 
time, an 0{1)- approximation to the optimal solution. 



Proof: Consider the associated hypergraph G = (V, E). By Lemma [4. 151 this hypergraph has 
the bounded growth property with Fk{t) = tO{k'^) (here 7(t) = 1). Therefore, the algorithm 
of Theorem 13.111 provides the required approximation. ■ 

Lemma 4.17 Given a set P of n points in the plane, and a parameter k, one can compute 
a set 1) of 0{k'^n log n) axis-parallel rectangles, such that for any axis-parallel rectangle b, if 
|b n P| < k, then there exists two rectangles bi, b2 G 2) such that (bi U b2) fl P = b fl P. 

Furthermore, consider the graph where two points of P are connected if they belong to the 
same rectangle in D. Then the number of edges in this graph is 0{nklogn). 
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Proof: Find a horizontal line i that splits P equally, and compute all the skyline rectangles 
that contain at most k points of P (that is, compute both the rectangle above and below the 
line). By Lemma [4.151 the number of such rectangles is 0{nk'^). Now, recursively compute 
the rectangle set for the points above £, and for the points below L Clearly, the number of 
rectangles generated is 0{k'^n\ogn), and let 2) denote the resulting set of rectangles. 

Now, consider any axis-parallel rectangle b such that |b fl P| <k. If it is does not intersect 
£ then by induction it has the desired property. Otherwise, if b intersects £, then it can be 
decomposed into two skyline rectangles, each one of them contains at most k points of P. 
By construction, for each of these rectangles there is a rectangle in D that contains exactly 
the same set of points. 

As for the second claim, we apply a similar argument. Consider an edge pq in this graph 
that arise because of a top skyline rectangle of L Furthermore, assume that p is higher than 
q and to its right. Clearly, there are at most k such edges emanating from p, as the skyline 
rectangle having p as its top right corner and having its left edge through q contains at most 
k points, and each such rectangle corresponds to a unique edge. As such, we get that the 
number of edges in the graph is E{n) = 0{nk) + 2T(n/2) = 0{nklogn). m 

Remark 4.18 A slightly more careful analysis shows that the number of rectangles in the 
set computed by Lemma [4. 171 that contain exactly k points is O {nk log n). This will not be 
needed for our analysis. 

Remark 4.19 Consider an instance G = (P, S) of PackPntsInRects. Let G' = (P, D) be 
a modified instance of PackPntsInRects where D is obtained from S by replacing each 
rectangle by two new rectangles whose union covers that same set of points. Lemma 14.171 
guarantees that this can be done such that |D| = O(n'^logn). One might be tempted to 
believe that we can plug G' into Theorem 13. Ill in order to get a bi-criteria approximation for 
G. Unfortunately, this does not work since (as the following example shows) the hypergraph 
does not have the bounded growth property for any meaningful growth function. 

Consider two parallel lines in the plane with 
positive slope. Place 6 (log n) points of P on each 
line (or close to the line) such that the points on 
(or close to) the top line all lie above and to the 
left of those on the bottom line. The remaining 
points of P will all lie below the points on the 
diagonals. Let T be the interval tree of P (using 
vertical split lines). Specifically, the remaining 
points of P will be placed such that each point 
on the top line lies in a different level of T, and all 
the points on the bottom line lie in the same node 
as the rightmost point on the top line. More 
specifically, the leftmost point on the top line will 
correspond to the root and the points in order 

from the left to right on the top line will correspond to continually walking down in the tree. 
(A @ in the figure represents a cluster of u points close together.) Now let X be the subset 
of P which consists of the two set of 0(log?7,) points on the diagonal lines. Consider the 
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intersection sub-hypergraph induced by X. Suppose that S has rectangle for every pair of 
points in P that can be obtained as the intersection of a rectangle with P. Then any pair 
of points from the top and bottom diagonals will correspond to a hyperedge in this induced 
sub-hypergraph. Therefore, Fi (log n) = Vt (log^ n) , and hence Theorem 13.111 only gives an 
0(£) approximation. 

Since (as the above remark demonstrates) we cannot directly apply Lemma 14.171 our 
approach will be more roundabout. We first show how to solve the independent set variant 
of our problem (i.e., unit capacities). Next, we slice the rectangles of the given instance 
with non-uniform capacities case into subrectangles with unit capacities, and plug it into the 
above algorithm to get a meaningful approximation. 

Lemma 4.20 Given an instance ofG = (V, 23) o/PackPntsInRects with unit capacities, 
one can compute a subset X C V, such that the total weight of X is r2(opt/log£) and each 
rectangle ofB contains at most 2 points of P, where n = |P|. 

Proof: We first use Lemma [273] to sparsify the given instance. We now have a set of P C V of 
t = G(£ log £) points, and an associated fractional solution, such that none of the constraints 
are violated. The value of the fractional solution on Gp is r2(opt), and as such we restrict 
our search for a solution to P. 

Furthermore, we can assume that the value assigned to each point of P by this fractional 
solution is exactly 1/M (we replicate a point i times if it is assigned value i/M), where 
M = 0(log£). Note, that none of the rectangles of !B contains more than M points of P. 
In particular, by Lemma [4.17[ one can build a set of rectangles 2) of size O(M^tlogt), such 
that every rectangle of S can be covered by the union of two rectangles of D; formally, for 
every b G S there exists bi, b2 G 2) such that bfl P = (bi U b2) fl P. We build a conflict graph 
G over P connecting two points if (i) they are both contained in a rectangle of S, and (ii) 
there is a rectangle of D that contains them both. By Lemma 14.171 this graph has at most 
0{Mt\ogt) = 0(£ log^ £) edges and t = e(£log£) vertices. 

We further add edges to G making a clique out of each group of duplicated points that 
arose from a single given point of P (this is needed since when duplicating the points we per- 
turbed them in order to maintain the implicit general position assumptions of Lemma I4.17[ 
and one needs to guarantee that at most one of these copies is picked to the independent 
set). Now for a point p G P with LP value Xp, the number of duplicated points is XpM. 
Hence the number of edges added for these cliques is 

pGP ^ '^ pGP peP 

and hence the number of edges in G overall is 0(£ log £). 

It is easy to verify that G has average degree 0(log^£), and the total weight of the 
vertices is G(optlog£), as such, by Turan's theorem, one can compute an independent set 
of vertices in this graph of weight r2(u7(P)/ (average degree -|- 1)) = r2(opt/log£). 

Now, it is easy to verify that any rectangle in 23 contains at most two points of this 
independent set. ■ 
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Theorem 4.21 Given an instance of {\/,'B) o/ PackPntsInRects (with arbitrary ca- 
pacities), one can compute in polynomial time a subset X C V that is an (0(log£),2)- 
approximation to the optimal solution. 

Proof: Compute a fractional solution to the given instance. Split each rectangle b with 
capacity #(b) into [#(b)/3] rectangles, each one containing at most value 4 from the frac- 
tional solution (this can be done by sweeping the rectangle from left to right, and splitting it 
whenever the fractional solution inside the current portion exceeds 3). Consider now a unit 
capacity instance on the same point set but with these new rectangles. We use Lemma 14.201 
in order to get an (0(log£), 2)-approximation for this new instance. 

We now show that this solution we obtained for the unit capacity instance is also an 
(0(log£), 2) -approximation to the original instance. First observe that the LP value on this 
new instance is r2(opt) (where opt is the LP value of the original instance) since scaling 
down the fractional solution to original instance by a factor of 4 would be a valid solution 
to the LP for the new instance (since these newly created rectangles each contained at most 
4 from the fractional solution), and hence the weight of the approximation is r2(opt/log£). 
Furthermore, we know every rectangle b G S contains at most 2 [7^(b)/3] < max(2, 7^(b)) 
points from this solution, since in the new instance each rectangle from S was replaced with 
[#(b)/3] rectangles each of which contains at most two points from the computed solution. 
(Note that the inequality holds since #(b) has integral value.) ■ 

5 Packing points into fat triangles 

In this section, we give a bi-criterion approximation for packing points into a set of a-fat 
triangles. More precisely, we consider the following problem. 

Problem 5.1 (PackPntsInFatTriangs.) Given a weighted set P of points and a set 7 of 
a-fat triangles in the plane such that each triangle A has a capacity #(A), find a maximum 
weight subset O C P, such that, for each A G T, the number of points of O contained in A 
is at most #(A). 

The approximation algorithm uses the following building blocks: 

(A) We prove that, for a given point set, there exists a small number of canonical sets 
such that for any fat triangle that covers at most k points, there exists a constant 
number of these canonical sets whose union covers exactly the same points. Showing 
this result is quite technical and requires non-trivial modifications of the work of 
Aronov et al. [AESlOj (in particular, their work does not imply this result). This is 
delegated to Section [5^ see Theorem 15.61 for the exact result. 

(B) An algorithm for approximating the unit capacity case. This follows by an algorithm 
similar to the one in Lemma r4.20[ see Lemma [5^ for details. Note that this uses the 
result from (A) to get the required approximation. 

(C) A partition scheme that shows that a fat triangle (with a measure defined over it) 
can be "partitioned" into 0{k) triangles such that any triangle in this partition has 
measure at most l/fc; see Lemma [5.31 
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Putting these components together yields the approximation algorithm; see Theorem 15.51 
for details. 

5.1 The unit capacity case 

Lemma 5.2 Given an instance G = (V, T) of PackPntsInFatTriangs with unit capaci- 
ties, one can compute a subset X C V such that the total weight of X is r2(opt/log^ £) and 
each triangle of 7 contains at most 9 points of P. 

Proof: We follow the proof of Lemma [4.201 We first use Lemma [2.31 to sparsify the given 
instance. We now have a set of P C V of t = G(£ log £) points and a corresponding fractional 
solution that is feasible. The value of the fractional solution on Gp is f2(opt), and as such 
we restrict our search for a solution to P. 

Furthermore, we can assume that the value assigned to each point of P by this fractional 
solution is exactly 1/M — we replicate a point i times if it is assigned value i/M — where 
M = 0(log£). Note that none of the triangles of 7 contains more than M points of P. In 
particular, by Theorem 15. 6[ one can construct a set Z of regions of size O (^M^t log t) such 
that, for every triangle of A G T, there exists a subset {zi, . . . , z^} (1 Z oi a.t most 9 regions 
(i.e. k < 9) such that P fl A = P n(u^^;^2;j). We build a conflict graph G over P connecting 
two points if (i) they are both contained in a triangle of T, and (ii) there is a set of Z that 
contains both of them. Since the number of sets in Z is O (^M^t log t) , and each such set has 
size at most M, it follows that the number of edges in the resulting graph G is O (M^t log^ t) , 
and the number of vertices is t = 0(£log£) vertices. As in the proof of Lemma [4.2U[ we 
also add edges between replicated points (since these edges do not affect our analysis, we 
ignore them for the sake of simplicity of exposition) . 

The graph G has average degree 0(M^ log t) = 0(log £), and the total weight of the 
vertices is B(optlog£). Therefore, by Turan's theorem, one can compute an independent 
set of vertices in this graph of weight f2(w(P)/(average degree + 1)) = r2(opt/ log^ £). 

Finally, it is easy to verify that any triangle in T contains at most 9 points of this 
independent set. ■ 

5.2 Covering a measure on a fat triangle 

At this point we would like to use Lemma [5. 21 in order to get a bi-criteria approximation for 
the case in which the capacities are arbitrary, as we did in Theorem 14.211 However, doing 
so directly proves more challenging for fat triangles than axis parallel rectangles. This is 
because our general procedure requires that, given an object x of a given type such that the 
total fractional value of the points in x is non-zero, we need to be able to decompose x into 
0(#(x)) smaller objects of the same type such that, for each smaller object, the fractional 
value of the points in the object is only a constant. This can easily be done for axis parallel 
rectangles by using vertical splitting lines (as was done in Theorem I4.2ip . but it is more 
challenging for fat triangles. However, the following lemma shows that such a decomposition 
is still possible for fat triangles. 
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Lemma 5.3 Let fi be a measure defined over the plane, and consider a fat triangle A. Then, 
for any integer k, one can cover A by at most 18k fat triangles, such that the measure of 
each of these triangles is at most fi{A)/k. 

Proof: To simplify the presentation we assume that /i(A) = 1. We recursively build a tree 
on A by partitioning the original triangle A into 4 similar triangles as shown in Figure [1] 
Each node of this tree corresponds to a triangle from this recursive construction. We stop 
the recursive partition for a node v as soon as the measure of the triangle A„ associated 
with it is at most 1/k. 

Once we have the tree, we select a set S of nodes of the tree as 
follows. We find the lowest node v in the tree such that the measure of 
its corresponding triangle is at least 1/k. We add the node v to S and 
we treat the measure inside the triangle corresponding to v as being 0. 
We repeat this process until the measure left uncovered is smaller than \/^ Figure 1: 
1/k, at which point we take the lowest node covering the remaining measure and we add it 
to S. We also add the root of the tree to S. Note that the set S contains at most k + 1 
nodes: since the initial measure is one, we added at most k nodes to S that are not the root. 

We also add all the nodes in this tree that are the least common 
ancestor (LCA) of a pair of nodes in S. Let S' be the resulting set 
of nodes. Now we can show that, for any tree T and any subset R 
of nodes of T, the set of all LCA nodes of all of the pairs of nodes 
in R has size at most |-R| — CJ. Therefore there are at most 1 5*1 — 1 
nodes in S' — S and thus the size of S' is at most 2|S'| < 2{k + 1). 

Let T be the set of triangles induced by the triangles corresponding \r Figure 2: 

to the children of the nodes of S' (i.e., every node of S" gives rise to four triangles). Consider 
the partition of the original triangle formed by T. It is easy to verify that every face in this 
arrangement has measure at most 1/k, and the number of faces of this arrangement, denoted 
by n, is at most 8A; (observe that since the triangles arise out of a recursive partition, a pair 
of such triangles is either disjoint or contained in each other). Furthermore, each face of this 
arrangement is either a triangle (which is a scaled and rotated copy of the original triangle), 
or the difference of two triangles where one contains the other (having this property is why 
we took the children of every LCA). We will refer to a face which is the difference of two 
triangles as an annulus face. Let n' be the number of annulus faces in this arrangement. 
Clearly, n' < \S'\ = 2k, as one can charge an annulus face to the node of S" that induced the 
hole in this face. 




^^Let u, V be the pair of nodes in R whose LCA has maximum depth. Let z be the LCA of u and v, and 
let R' ^ R — {u, v} U {z}. The pairs in R and the pairs in R' have the same set of LCA nodes and thus the 
number of LCA nodes of a set of size r satisfies the recurrence /(r) < /(?' — 1) + 1 and /(2) = 1. 
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We claim that an annulus face can be covered by the union of six 
translated and rotated copies of the original triangle. The easy case, 
is when the hole is a scaled and translated reflection of the outside 
triangle, see Figure |2l where three triangles are sufficient. The 
other case is when the hole is a scaled translated copy of the outer 
triangle, which by attaching to the hole three translated and rotated 

copies of the hole, gets reduced to the other case, see FigureO (The V/ rigure 3: 

property used here implicitly is that the outer triangle of the annulus can be partitioned into 
translated and rotated copies of the hole triangle, as the hole arises out of a recursive partition 
of the outer triangle.) 

As such, there are n — n' triangular faces in this arrangement and n' annulus faces. Thus, 
the original triangle can be covered by (ra — n') + Qn' = n + bn' < 8/c + 10/c < 18/;; translated 
and scaled copies of the original triangle that cover it completely, and no triangle in this 
collection has measure that exceeds \/k. ■ 

Remark 5.4 Given a weighted set of n points defining the measure inside the given fat 
triangle, the cover of Lemma [5.31 can be computed in O(nlogn) time. This requires using 
known techniques used in constructing compressed quadtrees, see |Harllj for details. 



5.3 The result 

Theorem 5.5 Given an instance of {\/, 7) of PackPntsInFatTriangs (with arbitrary 
capacities), one can compute, in polynomial time, a subset X C W that is (0(log^£) ,9)- 
approximation to the optimal solution. 



Proof: Compute a fractional solution to the given instance. For any triangle A in the 
plane, we denote by £(A) = ^pgvnA-^p ^^^ total mass of the fractional solution inside A. 
Next, we get a constant capacity instance out of (V, T) by replacing each triangle of T by 
a "few" triangles covering it, such that the total mass of the fractional solution inside each 
of these new triangles is at most c = 4 ■ 18 ■ 9. Formally, consider a triangle A G T, and 
let k = [#(A)/c]. If #(A) < c then there is nothing to do (as £(A) < #(A) < c), so we 
assume that #(A) > c. Applying the algorithmic version of Lemma [531 see Remark 15.41 we 
cover A with at most 
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triangles, where the total mass of the fractional solution inside each of them is at most 
£(A)A<#(A)/[#(A)/cl<c. 

Now, consider the generated instance with these new triangles, where each such triangle 
has capacity one. To this end, scale down the solution of the LP by a factor of c. Clearly, we 
now have a uniform capacity instance with an associated (valid) fractional solution having 
value r2(opt) (where opt is the optimal LP value for the original instance). Furthermore, any 
solution to this unit capacity instance, would correspond to a solution to the original instance 



(since we covered each original triangle with at most 18A; < 



#(A) 
2-9 



< #(A) new unit 



capacity triangles). Plugging this instance into Lemma [5^ yields the required approximation. 
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Specifically, every triangle A G T contains at most 9 [#(A)/(2 ■ 9)] < max(9, #(A)) points 
of the computed set of points. ■ 

5.4 Canonical decomposition for fat triangles 

In this section, we show that given a set P of ra points in the plane, and a parameter k, one can 
compute a set S of 0(/c^nlog n) regions, such that for any a-fat triangle A, if |A fl P| < /c, 
then there exists (at most) 9 regions in S whose union has the same intersection with P as 
A does. 

Our construction follows closely the argumentation of Aronov et al. |AES10j . However, 
our construction is (somewhat) different and (arguably) simpler since we are considering a 
"dual" problem to theirs. In particular, since modifying Aronov et al. |AES10] to get our 
result is not obvious, we present it here in detail. 

5.4.1 Initial setup 



To construct the set of regions, S, we will use an approach similar to that of Lemma [4. 171 As 
was observed in |AES10] . we can restrict our attention to axis aligned right triangles whose 
hypotenuse differs by no more than say one degree from —45°, as measured from the positive 
X-axis (i.e. it is near isosceles and faces to the right). In the following, let A be an arbitrary 
such triangle that contains at most k points. 

We first construct a two level interval tree on P, where the 
first level partitions the points based on their x-coordinate, and 
the second level based on their ^/-coordinate (and the splitting 
line for each node goes through the median point). Let v be 
the highest node in the first level of the interval tree whose 
corresponding split line, i, intersects A. Let A^ and A_5. 
denote the portion of A to the left or right of £, respectively. 
Also, let u be the highest node in the second level tree rooted at 
the left child of v whose corresponding split line, p, intersects 
A^, and let A\^ and A^ denote the portion of A^ above or 




Figure 4: Decomposing A. 



below p, respectively. (Note that we may assume that there 

exists split lines i and p that intersect A and A^_, respectively, 

since such regions that contain no points can be skipped). In the following, let p denote the 

point of intersection between i and p. See Figure HI 

We now construct sets of canonical regions, T_^, 7\^, and T^ such that for any choice of 
A there exists constant number of regions ri, . . . , r^ in 7^ U T^ U T^, such that A fl P = 
|Ji('^in P), and m < 9. 

We achieve this by showing that in each case (i.e., A\^, A^ and A^) the region r under 
consideration can be transformed into a polygonal region with a constant number of points 
P fo 

s. 

p is 0{nk^) 



of P (or orientations) defining its bounding edges, and whose intersection with P is the same 
as ro- We then show that the number of such regions needed for a particular choice 



and 



^■^For each point on the bounding edges, we wih need to specify whether it is inside or outside the canonical 
region. This can be encoded by a string of length c, where c is some constant bounding the number of defining 
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Figure 5: Handling A_>. 
(except maybe for the 



In the following, let P^, be the subset of points of P stored in the subtree rooted v, and 
let P„^^ be the set of points stored in the subtree rooted at u. 

5.4.2 Handling the right portion of the triangle (A^) 

For any A, we know that A^ will be a homothet of A, 
whose vertical edge lies on i, see Figure HI We now transform 
A_^ uniquely such that two points of P lie on its hypotenuse 
(or one point and the hypotenuse is at an angle of —46°) and 
one point of P lies on its bottom edge. Start by translating 
the hypotenuse towards the lower left corner of A^ (while 
clipping it to A_^) until it hits a point, pi. Next rotate the 
hypotenuse clockwise around pi until it hits a second point p2, 
or its orientation is —46° (as we rotate we modify its length so 
that one endpoint of the hypotenuse stays on i and the other 
on the base of A_^). Next translate the base of A^ straight 
upwards (while clipping it and the hypotenuse as to maintain 
a right triangle) until it hits a third point ps (which may be 
the same as the rightmost point out of pi and P2). Observe 
that the resulting region has the same intersection with P as A_ 
points on the boundary). See Figure |3 

We now bound the number of such resulting regions. As- 
sume that P2 lies to the right of pi (the other case is handled 
similarly). There are n„ = |Pi,| possible choices for p2. Now 
consider the horizontal line segment that connects i and p2. 
Rotate this segment clockwise around p2 (while increasing its 
length so that the other endpoint stays on i) until it hits pi, see 
Figure O We know that all the points we hit in this sweeping 
process lie in the computed region, and hence we can only have 
swept over k points before reaching pi (i.e. given p2 there are 
at most k choices for pi. If p2 does not exist we start with the 
triangle formed by i and a horizontal and —46° line through 
Pi). Now imagine translating the horizontal segment connecting p2 and i straight downward 
till we hit ps (while increasing its length so that its right endpoint stays on the hypotenuse 
defined by pi and P2). Again we know that all the points we hit in this sweeping process 
must be in our canonical region, and hence we can only have swept over k points before 
reaching P3 (i.e. given p2 and pi, there are at most k choices for P3). 

Hence there are 0{n^k'^) such canonical regions for the node v. Since P^- fl P„ =0 for 
any Vi, Vj at the same level in the top layer tree, summing across a given level gives 0{nk'^) 
canonical regions, where n = |P|. Thus, summing over all nodes in the top layer tree gives 
0{nk^ log n) such canonical regions overall. 




Figure 6: 



boundary points. Hence we can specify the inclusion or exclusion of the boundary points while only increasing 
the number of canonical regions by a factor of 2"^ = 0(1), and hence we will not need to worry about such 
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5.4.3 Handling the top left portion of the triangle ( A.^) 



Here we must consider two cases, based on the possible locations of p. If p ^ A^ (see 
Figure H]), we have a homothet of A whose bottom edge lies on p, and therefore we can 
argue as in the A_j. case, that this gives rise to 0{nu,vk'^) different canonical regions, where 
f^u,v = \^u,v\- Summing over all possible nodes u and v gives 0{nk'^\o^n) such canonical 
regions overall. 

Now suppose that p G A^, see Figure [71 In this case, we can 
extend Ar^ to get a homothet of A whose right side was cut off 
by / in order to get A^^ (A' in Figure El). Clearly, we have that 
Pu,v n Ai^ = Pu,v n A'. We can now generate a canonical region for 
A' in a similar fashion as the A_j. case, since it is just a homothet 
of A with its base lying on p, and then we can cut off the portion 
to the right of /. This would imply that we can generate 0{nu^vk'^) 
such canonical regions for the nodes u and v, and so overall there 
are O [nk"^ log^ n) such canonical regions. Figure 7: 




i 


I 


^;p\ 



Figure 8: 



5.4.4 Handling the bottom left portion of the triangle ( A^) 

Again we consider two cases, based on the possible locations of p. 
If p G A^ (see Figure (HI), then A^ is an axis parallel rectangle such p 
that one of its sides lies on £ (and another side lies on p). Hence by 
the proof of Lemma 14.171 in this case A^ gives rise to 0{k'^n\ogn) 
canonical regions overall. 

Now we consider (what is by far) the hardest case, when p ^ 
A^. In order to handle this case we will need to break up A^ as 
follows. Observe that A^ is a rectangular region whose upper right 
corner was cut off by the hypotenuse of A. First, we reduce A^ 
into a homothet of A, by removing rectangles r and r' from the left 
and bottom parts of A^, respectively (see Figure H}). This can be 
done since we already observed that by the proof of Lemma 14.171 
we can construct a set of O (n/c^ log^ n) canonical rectangles such Figure 9: 

that any rectangle (with a side on one of the split lines) has the same intersection with P as 
one of the canonical rectangles. For simplicity we continue to refer to the remaining part of 
A^ as just A^. 

We now break up A^ into three regions. Let /?"'" and /3~ denote 
the rays emanating from p at angles —140° and —130°, respectively 
(again, as measured clockwise from the positive a;-axis). These two 
lines split A^/ into three regions, which we will denote in their coun- 
terclockwise order as A^, A*/ and A"/ (see Figure ITOl) . Let a and 
h denote the intersection of /3+ with the hypotenuse and left edge of 
A^, respectively. Similarly, let c and d denote the intersection of 

(3~ with the hypotenuse and bottom edge of A /, respectively, 
rigure 10: 
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V' 



Figure 11: 
until it hits a third point, ps. 



Handling the top and bottom parts of A^ (i-e., 

A^ and A^). We now construct the canonical regions 
for A^. The construction is nearly identical to that for 
A^ and is included for the sake of completeness. The 
construction for A^ is omitted as it is symmetric to the 
Ay case. 

Start by translating the part of the boundary that in- 
tersects the hypotenuse of A towards the lower left corner 
of Ap (while clipping it to A^) until it hits a point, pi. 
Next rotate this edge counterclockwise around pi until it 
hits a second point p2, or its orientation is —44° (as we 
rotate we modify its length so that one endpoint stays on 
/S"*" and the other on the boundary Ay). Next translate 
the vertical edge of Ay to the right (while clipping it to A^ 

As for the number of such resulting regions, assume that p2 lies to the left of pi (the other 
case is handled similarly). There are nu,v = \Pu,v\ possible choices for p2. Now consider the 
vertical line segment that connects (3~^ and p2. Imagine rotating this segment counterclock- 
wise around p2 (while increasing its length so that the other endpoint stays on /J"*") until it 
hits pi. We know that all the points we hit in this sweeping process must be in our canonical 
region, and hence we can only have swept over k points before reaching pi (if p2 does not 
exist we start with the triangle formed by /3+ and a vertical and —44° line through pi). Now 
imagine translating the vertical segment connecting p2 and /?"'" to the left until we hit pa 
(while increasing its length so that its top endpoint stays on the line defined by pi and p2 
and its bottom endpoint on /?"''). Again we know that all the points we hit in this sweeping 
process must be in our canonical region, and hence we can only have swept over k points 
before reaching pa. Hence there are 0{nu,vk^) such canonical regions for a pair nodes u and 
V. Thus overall there are O (n/c^ log^ n) such canonical regions. 

^ P. 

Handling the middle part of A^ (i-e., 

A^). Let P°,, denote the subset of P„,„ 
that lies in between P~^ and f3~. Let the in- 
tersection of the hypotenuse of A with Ay 
be called the hypotenuse. 

Translate the hypotenuse towards p (while 
clipping it to A^) until it hits a point pi. 
Then rotate the hypotenuse clockwise around 
Pi until it hits a point p2, or it becomes verti- 
cal. Without loss of generality, assume that 
Pi lies to the right of p2. Let the horizontal 
(resp. vertical) line connecting pi and f3~^ 
(resp. /3~) be called h (resp. u). Translate 
h downwards (resp. u to the left), while en- 
larging it so that one endpoint stays on P~^ 
(resp. /3~), until either it hits the lowest 




Figure 12: 
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(resp. furthest to the left) point of A^ or a point outside of A^. Let this point be denoted 
P3 (resp. P4), and let h' (resp. z/') be the final translation of h (resp. u). See Figure [T2l 

Consider the region, r, bounded by the portion of h' to the left of p4, the portion of v' 
below ps , /S"*", /3~, and the line going through pi and p2 (this is the red shaded region in 
Figure [T2l) . First observe that if both ps and p4 lie outside of A^ then r will not cover all 
the points in Ay nP^^. Namely, the points lying in the rectangle defined by h\ u', and the 
vertical and horizontal edges of Ay might not be covered by r (see Figure W2\i . However, 
we already constructed a set of 0(/c^n log n) canonical rectangles, which we know contains 
two canonical rectangles that cover these points, and as such we do not have to worry about 
covering these points. Clearly, all the points of A°/ dP^^ either lie in this rectangle or in r. 
Next observe that there are no points of P^^ that lie in r that are not in A^ fl P^v- This 
follows from the easily proven fact (i.e. tedious but straightforward arguments) that since 
the hypotenuse was within one degree of —45°, that hn P^ lies to the right of b and u (1 f3~ 
lies above d. 

We now bound the number of canonical regions of type r. There are | -P^„ | possible choices 
for pi (which again we assume is to the right of P2). Now consider rotating h clockwise around 
Pi until we hit p2. We know from above that all the points we sweep past in this process 
must be contained in Ay r\ P^v ^^^ so given pi there are at most k possible choices for p2. 
Now consider translating h downward (resp. u to the left) until we hit pa (resp. P4). Again, 
from above we know that all the points we sweep over in this process must be contained in 
Ay n -P°^, and so given pi and p2, there are at most k possible choices for pa (resp. P4). 
Hence there are 0[k^ \P^^\) such canonical regions for a given pair of nodes u and v, and 
so overall there are O {k'^n log n) such canonical regions. 

5.4.5 Putting things together 

Summing the above bounds over all choices of the nodes u and v results overall in O [k^n log n) 
canonical regions. Furthermore, for any choice of A, we showed above that there exists a set 
of at most 9 of theses canonical regions whose (union of) intersections with P is the same as 
that of A. We thus get the following result. 

Theorem 5.6 Given a set P of n points in the plane, and parameters k and a > 0, one can 
compute a set S of 0{k^n\og ri) regions, such that for any a-fat triangle A, if |A n P| < k, 
then there exists (at most) 9 regions in S whose union has the same intersection with P as 
A does. 

6 PTAS for Unweighted Disks and Points 

In this section, we consider instances of the PackRegions problem in which the regions are 
disks with unit weights and all points have unit capacities. We now outline a PTAS for such 
instances based on the local search technique. The algorithm and proof are an extension of 
those of Chan and Har-Peled |CHlll[UH09] . and Mustafa and Ray [MRlOj . 

The algorithm. Since all of the regions have unit weight, we may assume that no region 
is completely contained in another. We say that a subset L of D is b-locally optimal if L 
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is a pointwise independent set and one cannot obtain a larger pointwise independent set by 
removing i < b regions of L and inserting i + 1 regions of D \ L. 

Our algorithm constructs a 6-locally optimal solution using local search, where b is some 
suitable constant. We start with L <f- 0. We consider each subset X C 2) \ L of size at most 
6+1: if X is a pointwise independent set and the set F C L of regions pointwise intersecting 
the objects of X has size at most \X\ — 1, we set L ^ (L\l^) UX. Every such swap increases 
the size of L by at least one, and as such it can happen at most n = |D| times. Therefore the 
running time is bounded by 0(?7,''+^6 |P|), since there are (j,"^) subsets X to consider and 
for each such subset X it takes 0{nb |P|) time to compute Y. 

Analysis. Let opt be the maximum pointwise independent set, and let L be the 6-locally 
optimal solution returned by our algorithm. If we can show that the pointwise intersection 
graph of opt U L is planar then the analysis in |CH11] will directly imply that |L| > (1 — 
0(lM))|opt|. 

We map the disks in opt and L to sets of points Qopt and Q\_ in IR^, respectively, and 
we map the points in P to a set of halfspaces Hp, by using the lifting of disks to planes and 
points to rays, and then dualizing the problem (see Section 1^^ . Mustafa and Ray prove that 
a range space defined by a set of points and halfspaces in IR^ has the locality condition, 
which is defined as follows. 

Definition 6.1 ( [MRIO] ) A range space R = (P, 2)) satisfies the locality condition if 

for any two disjoint subsets R,B C P, it is possible to construct a planar bipartite graph 
G = (-R, B, E) with all edges going between R and B such that for any D ^ T), if D n R ^ ^ 
and Do B y^ (/}, then there exist two vertices u ^ DCi R and v ^ Dn B such that {u, v) & E. 

Since opt and L are both pointwise independent sets, we know each point in P can 
intersect at most one disk from opt and at most one disk from L. Hence each half space in 
Hp can contain at most one point from Qopt and at most one point from Qi. Since points 
and halfspaces in IR^ have the locality condition, setting R = L and B = opt immediately 
implies that there is a planar graph on the vertex set L U opt such that any vertex from 
L and any vertex from opt that are in the same halfspace are adjacent. In particular, the 
intersection graph is planar. 

Theorem 6.2 Given a set of n unweighted disks and a set of m points in the plane (with 
unit capacities) , any b-locally optimal pointwise independent set has size > {l—0{l/\/b))opt, 
where opt is the size of the maximum pointwise independent set of the disks. In particular, 
one can compute an independent set of size > (1 — £)opt, in time mrP^^^'^ ^ 

Corollary 6.3 There is a PTAS for instances o/ PackHalfspaces in which each halfspace 
has unit weight, and each point has unit capacity. 

Corollary 6.4 There is a PTAS for instances of PackRegions in which each region is a 
unit-weight disk, and each point has unit capacity. 

Corollary 6.5 There is a PTAS for instances of PackPoints in which each region is a 
unit-capacity disk, and each point has unit weight. 
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7 Hardness of approximation 

7.1 Packing same size fat triangles into points 

Here we show that PackRegions (Problem ll.2p does not have a PTAS, even if the regions 
have unit weight and their union complexity is linear. We show that the problem is APX- 
hard using a reduction from the maximum bounded 3-dimensional matching problem. Since 
maximum bounded 3-dimensional matching is APX-complete |Kan91] . this will imply the 
claim (unless P = NP). 

Theorem 7.1 Unless P = NP there is no PTAS for PackRegions (Problem lTT^) even if 
the regions are unweighted, in the plane, and have linear union complexity. In particular, this 
holds if the regions are fat triangles of similar size. (See Corollary \4.!^ (EP for the matching 
approximation algorithm.) 

Proof: Let T(lAxBxChe the input triples for an instance of maximum bounded 3- 
dimensional matching, where A, B, and C are disjoint subsets of some ground set X (for 
simplicity we assume X = AU B U C). For each element x G X we make a representative 
point Vx and place it arbitrarily on the unit circle in the plane and give it unit capacity. Let 
Va, Vb, and Vc be the sets of representatives for A, B, and C (respectively). A triple in T 
thus corresponds to a triangle with one vertex in each of Va, Vb, and Vc- Clearly, finding 
a maximum packing of these triangles into these points is an instance of PackRegions. 
Moreover, a maximum packing here corresponds to a maximum set of triangles (triples) 
such that each point (element of X) is covered by at most one triangle. Therefore a PTAS 
for this problem translates to a PTAS for the maximum bounded 3-dimensional matching 
problem. (Note that this does not imply that there is no PTAS for other specific types of 
regions.) 

Now we show that we can make the triangles fat and of similar size, 
and hence there is no PTAS even in the case of linear union complexity. 
Let the range of a set of representative points be the angle around the 
circle between the farthest two points of the set, and let the center of a set 
be the midpoint on the circle between the farthest two points of the set. 
Instead of placing the points arbitrarily, we will place the points so that 
the range of each of Va, Vg, or V^ is less than five degrees. Moreover, we 
place the points so that the centers of Va, Vb, and Vc are 120 degrees apart. In this case the 
triangles all have roughly the same size and are nearly equilateral. It is known that such a 
set of triangles has linear union complexity |MMP"'"94j . Hence, by the above reduction, even 
in this case where the regions are restricted to have linear union complexity (and even more 
specifically when they are restricted to be fat triangles of roughly the same size) , we cannot 
get a PTAS. ■ 

7.2 Packing points into fat triangles 

Lemma 7.2 There is an approximation-preserving reduction from the INDEPENDENT Set 
problem in general graphs to the PackPoints problem. In particular, for instances of the 
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problem PackPoints in which the regions are fat triangles with unit capacities and the 
points are unweighted, no approximation better than Q{n^~'^) is possible in polynomial time, 
for any constant 5 > 0, unless P = NP. 

Proof: Consider an instance of the INDEPENDENT Set problem, namely a graph G = {V, E). 
Let n = \V\. Place n distinct points on the unit circle (arbitrarily) and map every vertex of 
\^ to a unique point of the resulting set of points P. For every edge uv G E, consider the 
segment PuPv, where p„ and p^ are the points corresponding to u and v in P. We construct 
a fat triangle containing p^p^ by connecting p„, p^, and a third vertex in the interior of the 
unit disk; this can always be done so as to achieve roughly 2-fatness. We add this triangle 
to our set of regions 2), and assign it capacity one. 

Clearly, solving the resulting instance (P, D) of PackPoints is equivalent to solving the 
Independent Set problem for G. The claim now follows from the hardness results known 
for the Independent Set problem |Has99j . ■ 

8 Conclusions 

In this paper, we presented a general framework for approximating geometric packing prob- 
lems with non-uniform constraints. We then applied this framework in a systematic fashion 
to get improved algorithms for specific instances of this problem, many of which required 
additional non-trivial ideas. There are several special cases of this problem for which we 
currently do not know any useful approximation; for example, the special case of packing 
axis-parallel boxes into points, in which the boxes are in four dimensions is still wide open. 
Making some progress on these special cases is an interesting direction for future work. 
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